Quantum radiation pressure on a moving mirror at finite temperature 
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We compute the radiation pressure force on a moving mirror, in tlie nonrelativistic approxima- 
tion, assuming the field to be at temperature T. At high temperature, the force has a dissipative 
component proportional to the mirror velocity, which results from Doppler shift of the reflected ther- 
mal photons. In the case of a scalar field, the force has also a dispersive component associated to a 
mass correction Am. In the electromagnetic case, the separate contributions to the mass correction 
from the two polarizations cancel. We also derive explicit results in the low temperature regime, 
and present numerical results for the general case. As an application, we compute the dissipation 
and decoherence rates for a mirror in a harmonic potential well. 
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In his seminal paper published in 1948 [1], Casimir computed the attractive force between two neutral perfectly 
O , conducting plates due to vacuum fluctuations of the electromagnetic field. The Casimir force itself is a fluctuating 
quantity [2], and from the general argument related to the fluctuation-dissipation theorem [3] one may expect that 
dissipation occurs in the case of moving boundaries. The energy dissipated from macroscopic moving bodies yields 
for the creation of real particles (photons in the case of the electromagnetic field) [4] . Hence the vacuum radiation 
I pressure on moving boundaries has a dissipative component that plays the role of a radiation reaction force. 
J> , This effect takes place even in the case of a single moving plate, as shown by Fulling and Davies [5]. They treated 
' exactly the problem of a massless scalar field in 1-1-1 dimensions in the presence of a plate moving in a prescribed 
, arbitrary way. However, since they employed a method based on conformal transformations, their results could not 

■ be generalized to higher dimensions. In order to address the case of 3-1-1 dimensions in the non-relativistic regime, a 
\ convenient perturbative method was proposed by Ford and Vilenkin [6]. Their approach is based on the assumption 

■ that the field modification induced by the motion of the plate is a small perturbation, which is computed up to first 
order on the displacement of the plate. They considered a massless scalar field, in either 1-1-1 or 3-1-1 dimensions. 
In the former case, the dissipative force is proportional to the third time derivative of the plate's displacement, and 

'"pi ' corresponds to the non-relativistic limit of Fulling and Davies' result. For 3-1-1 dimensions the force on a plane mirror 
' , moving along the normal direction is proportional to the fifth time derivative of the displacement. This is also the 
case when the electromagnetic field is considered [7] , although the proportionality factor is not simply twice the value 
found for the scalar case, as would be guessed by crude analogy with the static Casimir effect. Higher order derivatives 

■ appear when considering moving mirrors of finite extent [8] . 
A dissipative force proportional to the velocity of the mirror (like a viscous force) would clearly violate the Lorentz 

invariance of the vacuum field. For a thermal field, on the other hand, this requirement docs not hold, and the thermal 
?H ' contribution to the dissipative force turns out to be proportional to the velocity in the case of 1+1 dimensions [9]. 
The effect of thermal photons is larger than the contribution of vacuum fluctuations to the force for temperatures 
larger than hujo/kB (where ks is the Boltzmann constant, and ujq a typical mechanical frequency). This corresponds 
to temperatures in the mK range for frequencies in the MHz range. This clearly shows the importance of temperature 
in the dynamical Casimir effect, which would probably provide the dominant contribution in any attempt to measure 
the force. Thermal effects on the generation of photons in a cavity with moving mirrors have also been considered [10]. 

In this paper, we analyze the thermal contributions to the radiation pressure force in 3-1-1 dimensions, for both 
scalar and electromagnetic fields. We take a perfectly-reflecting plane mirror moving along the normal direction, in 
the non-relativistic regime. Our approach allows us to identify and distinguish between the fleld modes contributing 
to the dissipative component of the force from those contributing to its dispersive component. We derive analytical 
results in the low and high temperature limits, and also compute numerically the force in the general case. The paper 
is organized as follows: in the next section we take a massless scalar field under Dirichlet boundary condition. In 
Sec. HI, we consider the electromagnetic field. The results of this section are then applied, in Sec. IV, to the analysis 
of dissipation and decoherence of a mirror in a potential well. Section V presents an interpretation of the results in 
the high-temperature limit and some concluding remarks. 
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II. MASSLESS SCALAR FIELD 



We choose Cartesian axis such that the plane of the mirror is parallel to the OXy plane. The mirror is displaced 
along the OZ direction in a prescribed, non-relativistic way. Hence the field </> satisfies the wave equation and the 
Dirichlet boundary condition; 

52^ = with (j){x,y,5q{t),t) = , (1) 

where 5q(t) denotes the position of the mirror at time t. We assume that 5q(t) is small when compared with some 
characteristic field wavelength. We follow the perturbative approach of Ford and Vilenkin [6] and write the field as 

^ = ^°+6cl>, (2) 

where (fP is the solution of the corresponding static problem: 

d^cjP = with / (a;, y, 0, i) = , (3) 

and 6(j) is a small motion-induced perturbation. By taking the Taylor expansion around 2: = up to first order in 6q, 
we derive the following boundary condition: 

S<i>{x, y, 0, t) = -Sq{t) d,clP{x, y, 0, t) + 0{5q^) . (4) 

We solve for 5(j} in terms of cfP in the Fourier representation, defined as (using capital letters for Fourier transforms) 

$(^,k||,a;) = j dt j d^ry e^'^'e"*'""-''" 0(r|| + zi,t). (5) 



We find 



where (we take c = 1) 



<5$(2,k||,a;) = -e^'=l-l j ^<50(u; - u;i„)a,$°(0, ky , a;;^). 



(6) 



K = [(w + ief - with e ^ 0+ 

is defined, for a given value of fcy, as a function of uj with a branch cut along the real axis between — fcy and /cy, so 



that K is positive for uj > k\\, negative for uj < —fey, and equal to *y /^y — otherwise. Then, when corresponding 

to a propagating field, propagates outwards from the region around the moving mirror; otherwise it corresponds 

to an evanescent wave. According to (6), the scattering by the moving plate generates frequency modulation: for a 
given mechanical frequency the input field Fourier component at frequency Win is scattered into a new frequency 
w = Win + t^o- Due to translational symmetry along the OXy plane, all scattered components have the same ky. If 
I Win + wol < fc||, the scattered wave is evanescent. 

We write the Fourier representation of the unperturbed field in the half-space z > in terms of the bosonic operators 
a(k): 



$0(2,ky,w)=iji^^^^sin(|K|^)[a(|K|z + ky)e(w)+a(|K|z-ky)te(-w)]e(K2), (7) 



where © denotes the step function. This equation shows explicitly the association between positive (negative) fre- 
quencies and annihilation (creation) operators. Moreover, the normal mode decomposition includes only propagating 
waves (since evanescent waves do not satisfy the required boundary condition), hence the factor G(k^). A similar 
representation may be written for the field in the half-space 2; < in terms of independent bosonic operators 6(k) 
(and their Hermitian conjugates). 
We compute the radiation pressure force from the stress tensor component 

T,, = i [{d.^f + {dy^f - {d.ct^f - iMf] , (8) 
taken at the surface of the moving mirror. Up to first order in 6q, the force is given by 
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Sf{t)= J dxdy {dT,,{x,y,0+,t) - dT,,{x,y,0- ,t)) 
where ({...,...} denoting the anti-commutator) 

5T^^{x,y,z,t) = {d^(lP{x,y,z,t),d^5(l){x,y,z,t)} 



(9) 



(10) 



is the motion-induced modification of the stress tensor. Its Fourier representation may be computed from Eqs. (6) 
and (10). When taking the average over a given field state we find (with e — > 0+) 



where 



(5T..(0+ k|„.)) =y'^Jj0 [(c.-c.,+^e)^-(k„-k,„)^]^/^ 
X J 5(3(a; - wi - W2) cr(ki||,a;i;k|| - ki||,W2) , 

a(l;2) = ({9,$0(0+,ki||,wi),9,$"(0+,k2||,c^2)}) 



(11) 



(12) 



is the correlation function of the unperturbed field taken at the OXy plane. For a thermal field, we find, using the 
normal mode decomposition as given by (7): 



a(l;2) = len^h^Lvf - 6(0)? - fcy?) <5(a;i + Wa)^^'^^! || + kay) [1 + 2n(|a;i| 



(13) 



where 



n{\uj\) = (e 



ft|w|/(feBT) 



1)" 



is the average photon number at temperature T. 

Wc now analyze in detail the expression in the r.-h.-s. ofEq. (13). The factor l-|-2n(|ciJi|) originates from the general 
relation between the thermal averages of anti-commutators and commutators, as given by the fluctuation-dissipation 
theorem [11]. For the field operators themselves, the commutator is a c- number, hence the temperature dependence 
comes solely from this factor. The factor (5'^^)(ki[| +k2||) is a signature of a homogeneous field state: it means that the 
correlation function for two given points on the surface of the plate depends only on the relative position between the 
points. When replaced into (11), it yields for an uniform pressure over the surface of the plate. Likewise, the factor 
5{u)i + CO2) in Eq. (13) is a signature of a stationary field state. When written in the time domain, it corresponds to 
a correlation function depending only on the time difference, not on the individual times themselves. When replaced 
in (11), this factor singles out the mechanical Fourier component at the same frequency lu appearing in the argument 
of (7^2(0+,k||,a;)). To compute the force from (9), we also need the motion-induced stress correction at 2 = 0~, 
which is computed from (6) and the normal mode decomposition of the field in the half-space 2; < in analogy with 
the derivation of (11). Its contribution simply doubles the value of the net force, which we write as (we employ the 
superscript D to denote the results for the scalar field obeying Dirichlet boundary conditions) 



5F°(a;)=x°(a;) SQico). 



(14) 



We replace Eq. (13) into Eq. (11) and integrate over the OXy plane to derive the susceptibility function X^if^) (^ 
is the area of the mirror): 



ihA f /■I'^'nl 

X°(w) = J dlOin J k 



1/2 



2 + '^(1'^" 



(15) 



where we have replaced the variable of integration wa in (H) by Win, since it corresponds to the frequency of the 
unperturbed field $0 [or input frequency, sec Eq. (6)]. 

The imaginary part of the susceptibility provides, according to Eq. (14), a force component in quadrature with the 
displacement. If Imx^ > 0, this force component is in opposition of phase with respect to the velocity of the mirror, 
and hence dissipates its mechanical energy. On the other hand, the real part provides the dispersive force component, 
which is in quadrature with the velocity, and does not engender any energy exchange when averaging over a sufficiently 
long time interval. According to Eq. (15), Rex° results from contributions of input modes satisfying [a; -|- Win| < k\\, 
i.e. (propagating) modes that generate evanescent waves when scattered by the mechanical Fourier component w. In 
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Fig. 1, we represent the region for the integration in the r.-h.-s. of (15), corresponding to the condition |wi„| > fcy > 
(all input modes are propagating waves), in the plane win — fcy- It is divided in four subsets, labeled Ri to R4 [7]. In 
this diagram, the scattering by the mechanical frequency lu corresponds to an horizontal displacement, by an amount 
equal to uj, from the point of coordinates (a;i„,A:||) representing a given input field mode (for the sake of clarity we 
assume w > in the diagram). The contribution to the dispersive component of the force comes from region i?4, the 
region that occupies the evanescent sector when shifted by uj, whereas Ri to R3 contribute to dissipation. 




-2(0 -co (0 2co 



FIG. 1. Diagram for the evaluation of the susceptibility, given by the integrals in the r.-h.-s. of (15) (scalar or electromagnetic 
TE-polarized modes) and (23) (electromagnetic TM-polarized modes). Field modes propagating along different directions 
corresponding to the same values of fc|| = |k||| and frequency are represented by a single point in the tJin — fcy plane of 
integration. Evanescent input modes (|cJin| < fcy, grey region) are excluded. Regions Ri to R3 yield the dissipative component 
of the force, whereas region R4,, corresponding to input propagating waves that are scattered into evanescent waves, provides 
the dispersive component. At zero temperature, the contributions from Ri and R2 cancel, because the integrand in these 
regions is anti-symmetric with respect to reflection around the axis at ojin — ~i^/2, which is indicated by a dashed line. 



There are two distinct terms in the r.-h.-s. of (15), both contributing to the dispersive and dissipative components: 
one proportional to ndwinl), corresponding to thermal fluctuations, and one independent of temperature, containing 
the effect of vacuum fluctuations (with n replaced by 1/2). Accordingly, we write the susceptibility as 

D D I D 

X = Xt + Xvac- 

At zero temperature, we have by definition Xt=o ~ ^ ^^'^ ~ X?ac- It is then particularly useful to consider 
the reflection around the axis t^in = —uj/2 (indicated by a dashed line in Fig. 1), which is implemented by the 
transformation u-m ^ -~<-^ — while keeping fc|| unchanged. The contributions from points in region R2 cancel 
exactly those from their conjugates in region i?i in the integral in Eq. (15). As a consequence, the single contribution 
to dissipation at zero temperature comes from R^, the only bounded region in the diagram: it corresponds to negative- 
frequency input modes that are scattered into positive-frequency propagating modes. We evaluate the resulting 
integral to find 

hA 

in agreement with Ref. [6]. Thus, the dissipative force exerted by the vacuum field is caused by the motion-induced 
mixture between positive and negative field frequencies. The discussion following Eq. (7) indicates that this mixture 
is a signature of a Bogoliubov motion-induced transformation of creation into annihilation operators (and vice- versa) , 
which is clearly connected to the emission of particles [12] [13]. In fact, the dissipative force in vacuum plays the 
role of a quantum radiation reaction force, dissipating the mechanical energy at exactly the rate required, by energy 
conservation, for the photon emission effect. 

As for the dispersive component of the force, on the other hand, the integral runs over the unbounded region _R4, 
and the vacuum contribution diverges. After regularization, the dispersive force leads to renormalization of the mass 
of the mirror [6], an effect analyzed in detail for a dielectric interface in Refs. [14] and [15] and for a dispersive mirror 
in Ref. [16]. 
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We now analyze the thermal contribution to the force, represented by %?• contrast with the vacuum force, 
the thermal dispersive force, as given by the integral over region i?4, is finite, because the average photon number 
decreases exponentially to zero at high frequencies. In Appendix A, we show that Rexri'^) < for any uj (on the 
other hand, the imaginary part must be positive so as to provide dissipation). The ratio x^/I^Xvac ^ function of a 
single parameter, ksT / [Tllo). Numerical results for |Rexy l/Imx^^^j. are shown in Fig. 2 (dashed line). The dominant 
contribution comes from input frequencies satisfying jcjinl ~ kT/h. When ksT <^ huj, this condition implies \uJin\ <C clJ, 
and hence the dominant contribution to Xt comes from the close neighborhood of t^in = in Fig. 1. However, region 
R4 is separated from this neighborhood (see fig. 1), and the larger contribution comes from nearly grazing field modes 
with frequencies close to ujin — — The corresponding average photon number is n ~ exp(— ?i|w|/2fcBT), hence 
Rex^ is exponentially small in this limit. In Appendix A, we derive 

^ n . ^ 1 {ksTf , n / h\uj\ \ , , 

Rex?.o(-) = -^4^^..^exp (-^j . (17) 

Fig. 2 suggests that |Rex^| grows according to a power law when kBT/[fiLo) ^ 1. Neglecting terms of the order 
of {Tioj /ksTY , we derive in Appendix A the following result for the dispersive thermal susceptibility in the high- 
temperature limit: 

Rex?_oo(-)^-^^^^c.^ (18) 
where ^(3) ~ 1.2, C denoting the Ricmann zcta function [17]. 
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FIG. 2. Thermal susceptibility for the scalar field. The solid line represents Ihixt, whereas the dashed line represents —Rex?, 
both divided by ImXvac- In the insert, we plot Imx? divided by its high-temperature asymptotic value, showing the cross-over 
between the low and high temperature regimes. 

The thermal dissipative susceptibility Imxj^ is computed from Eq. (15) in a similar way. We define [with t = 

kBT/{nuj)] 



D _ Imxg(^) 



(19) 



In Appendix B, we derive the results G°(t) = 2At:'^t^ for r > 1, and G°(t) = I&'k'^t^ for t < 1. We plot G°(t) 
in Fig. 2 (solid line), showing that the thermal contribution to dissipation becomes larger than the vacuum effect for 
UbT / {Two) ~ 0.2. However, the deviation from the high-temperature behavior is not visible in this plot. In the insert of 
Fig. 2, we plot the ratio Imx?('-^)/IiiiXT-+oo('-^)i showing the smooth cross-over between the low and high temperature 
regimes. A similar cross-over occurs when considering the dissipative susceptibility for the electromagnetic field, as 
discussed in the next section. 
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III. ELECTROMAGNETIC FIELD 



We consider the following boundary conditions for the electric and magnetic fields E' and B' measured in the 
instantaneously co-moving Lorentz frame S' : 



We write 



z X E'(a;', y',z' = 0) = z- B'{x', y', z' = 0) = 0. 



B{z, k|| , c^) = E(™) k|| , to) + E(™) k|| , u;) , 



(20) 



where E^"""^) is perpendicular to the plane defined by the vectors z and ky, whereas E^'^'^) is parallel to this plane. 
When considering the scattering by the moving plane mirror, these two polarizations are not mixed, and may be 
mapped into two scalar-field boundary value problems. For TE polarization, Eq. (20) yields a Dirichlct boundary 
condition for the vector potential in the laboratory frame identical to Eq. (1). The contribution of TE-polarized field 
modes coincides with the results found for the scalar field in the previous section: x^™^ = X^- 

In order to compute the contribution from TM-polarized modes, we follow the approach of Ref. [7] and define a 
new vector potential: 



E(™)=Vx^; B(™) = 



dA 
'dt' 



Prom Eq. (20), we derive a Neumann boundary condition for A. in the co-moving frame [18]: 



dA 

dz' 



j{x',y',z' = 0) = 0. 



(21) 



(22) 



From this point we proceed as in the previous section. After a lengthy calculation, we find for the contribution of 
TM-polarized modes 



X 



(TM) 



= J dujin J dk\\ k 



kf. - UJinioJ + Wi, 



(w -h Win + ie)^ - ku 



V2 



-k^ 



+ n{\ujin\) 



(23) 



The integration region for the evaluation of 

^(TM) ig 

divided as shown in Fig. 1, with R4 providing its real part, and 
Ri to i?3 its imaginary part as in the scalar case. The vacuum contribution was already discussed in detail in Ref. [7]. 
For the imaginary part, the contributions from regions Ri and R2 cancel, and R3 yields a contribution larger than 
the TE result, so that the total dissipative susceptibility for the electromagnetic case is not simply twice the result of 
Eq. (16) for the scalar field: 



Imx^,':? = Imxl™)+Imx(™) 



hA 

30^' 



(24) 



The TM contribution to the real part, on the other hand, cancels the (inertial) term from TE modes, but a 
divergent term remains [7]. 

A similar cancelation takes place when considering the thermal contribution in the high-temperature limit. Starting 
from Eq. (23), we show in Appendix A that Rex^'^^ {u>) is positive for any u> (whereas the TE contribution is negative), 
and derive 



Rex?_!l(-) 



C(3) {kBTf 
2n 



(25) 



Thus, the TM contribution cancels the TE contribution as given by (18) to leading order of kBT/{hw), and the 
resulting electromagnetic dispersive susceptibility is smaller than the scalar susceptibility by a factor of the order 
of {Tiu/kBTY . For finite values of kBT/{hLo), Rexj?'^^('^) is larger than |Rex7^^^(w)|, so that Rex™(w) is always 
positive. In Fig. 3, we plot Rexj^ divided by — RexT^oo('^) ^ ^ function of kBT/{Jiui), showing the {kBT/Jkj)~^ 
behavior at high temperatures. 

In the low temperature limit, on the other hand, both contributions become exponentially small. In Appendix A, 
we find 
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Rex?^o(-) - ^ksTA^'c^p (^1^) ' ^^^^ 

The leading contribution to both TE and TM terms originate from field modes propagating close to grazing directions, 
and with frequencies close to —uj/2. As discussed in Sec. II, the resulting susceptibility is exponentially small because 
it is proportional to the average photon number at frequency |w/2|. As compared to TE field modes [see (17)], TM 

modes provide a larger contribution, by a factor of the order of {Tilj /ksTY , so that Rexf^o('^) ~ ^^Xr-^oi'^)- This 
is in line with the preference for TM polarization when considering propagation along a direction parallel to the plane 
of the mirror, since it matches the boundary conditions in the static case. 
For the imaginary part, we define, in analogy with (19), 



As discussed in Appendix B, TM and TE modes provide identical contributions in the high temperature limit, hence 
Imx™ is simply twice the value for the scalar field in this limit: G^^{t ^ 1) 487r^r^. The TM contribution in 
the low temperature limit is also analyzed in Appendix B. We find G^^{t <C 1) ~ 647r'*r^. Numerical results are 
presented in the next section, where we analyze dissipation and decoherence in a harmonic potential well. 




kgT /(ti(o) 



FIG. 3. Dispersive thermal susceptibility for the electromagnetic field. We plot Rex™/|RexT^cx) I a-s a function of ksT /{hui). 



IV. DISSIPATION AND DECOHERENCE 



In this section, we consider the effect of the quantum radiation pressure force on the motion of the mirror. We 
start with a classical description of the position of the mirror, and analyze the damping of the mirror's oscillation in a 
harmonic potential well (frequency wo)- Later in this section we also consider the quantum dynamics of the mirror, in 
order to derive the decoherence rate induced by radiation pressure. The limiting cases of zero and high temperatures 
were considered in Refs. [19] and [20], respectively. The results of the previous section allow us to take arbitrary 
values of temperature. 

The equation of motion reads 

m5q{t) = -mwl5q{t) + j dt'x{t')Sq{t - t'), (28) 

where x{t) denotes the inverse Fourier transform of xi'^)-! with x{^) = X^{^) ^or the scalar field model and x(^) — 
X^^{u!) for the electromagnetic case. Its solutions arc of the form 6q{t) = e^^^^*6qa, where the complex constant is 
found by replacing this function into (28): 

-mD,'^ = -mujl + xi^)- (29) 
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We include the dispersive effect of the vacuum field through a renormalization of the bare mass and frequency 
of oscillation, so that m and ljq in (28) and (29) are already renormalized and x iii (29) is replaced by Rcxt + 
i(ImXvac + Ihixt)- This allows us to assume that the effect generated by the radiation pressure force is a small 
perturbation of the free oscillations at frequency loq when computing the roots of Eq. (29). Hence they are of the 
form O = ±(a;o + SijJr) — «r, with F, \5ujt\ "C wq- F is the damping rate, and Slot is the frequency shift induced by 
the coupling with thermal photons. Using that Rexr (Imx) is an even (odd) function of f2, we find from (29) 

RexT(^o) ,„„^ 

SujT = (30) 

2mLOo 

and 

F=^^^^. (31) 

2ma)o 

The case of a mass correction Am provides a trivial application of (30) (this is the case for the 3D scalar field for 
frequencies much smaller than fc^T/Tl, but not for the electromagnetic field): when the dispersive susceptibility is of 
the form RexT(w) = Am w^, (30) yields Slot/loq = —Am/ {2m), which is just the first order expansion of the exact 
result {loo + Slot)/loo = \/m/{m + Am). 

As discussed in Sec. II, lmx{too) is always positive since energy is taken from (and not given to) the mirror; hence 
F, as given by (31), is positive, which corresponds to exponential decay as required. 

Of particular interest is the application of (31) to the electromagnetic case, which we analyze in the rest of this 
section. At zero temperature, Eqs. (24) and (31) yield (here and in the next section we re-introduce the speed of light 
c) 



307r2 2mc2 

Hence Fvac is proportional to the ratio between the zero point and the rest mass energies of the mirror. Despite of the 
geometrical factor AlOq/c^, representing the squared ratio between the transverse size of the mirror and the typical 
vacuum field wavelength (for frequencies in the GHz range c/loq is in the centimeter range), we have Fyac ^ ^o, as 
required for consistency of the derivation. 

In the high-temperature limit, Eq. (31) yields 

r^_^^M:A(3z)!, (33) 

15 mc^ ' 

in agreement with Ref. [20]. Except for a numerical factor, (33) differs from (32) by the replacement of the zero point 
energy %loo/2 by ksT and the frequency scale loq by kBT/U. In Fig. 4, we plot 

^ Imx^^ = 1 , 1 rEM fkBT\ 
Fvac Imx^a'^ 12 \hoooJ 

as a function of ksT/ (Hlo) (solid line), showing the T** behavior corresponding to the high-temperature approximation 
as given by Eq. (33), which is only 4% larger than the numerical exact value at kBT/{hLOo) = 1. 

In order to analyze the effect of decoherence, we now consider the quantum dynamics of the mirror in the 
potential well. We assume that the mirror state is initially a coherent superposition of two coherent states: 
IV-) = (|ao) + I - ao)) /2, with 

1 AZ 



where AZ represents the distance between the wavepacket components at t = 0, and AZq = ^Jh/2m(jjQ is the 
uncertainty of position of the ground state. 

The free evolution preserves the form (and the coherence) of the initial state, with ao replaced by a{t) = 
ao exp(— iwoO- The very weak radiation pressure coupling with the field introduces two new time scales, both as- 
sumed to be much larger than the free evolution time scale 27r/wo. The time needed to reach thermal equilibrium 
is the largest time scale. At zero temperature, it corresponds to l/(2r), where V is the amplitude damping rate 
discussed before in the classical theory. In fact, in this particular case, the equilibrium state is of course the ground 
state of the harmonic oscillator. During the decay process, the energy quantity (2aQ — 1/2)?ujOq is dissipated by the 
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dynamical Casimir effect, with the emission of about 2a^ pairs of photons, since the energy contained in each pair is 
hujQ [13]- [15] [21]. 




FIG. 4. Damping and decoherence rates for a mirror in a harmonic potential well. We plot F/Fvac (solid line) and F'*'"^/F' 



dec 

vac 



(dashed line) as functions of ksT / (huo) ■ Results for the damping and decoherence rates at zero temperature Fvac and F^^J; are 
discussed in the text; note that Fvac "C F^ac- 



The second time scale corresponds to the process of decoherence. The state decays into the incoherent statistical 
mixture p ^ {\a) {a\ + \ — a) (— a|) /2 much before reaching thermal equilibrium. The decoherence time is much shorter 
than the damping time because a single photon pair already contains 'which way' information sufficient to destroy 
the possibility of interference between the two components. Since the time scale corresponding to the emission of 2aQ 
pairs is l/(2r), the decoherence time at zero temperature is l/(4aor) = 4(AZo/AZ)^(l/r) <C l/F. The decoherence 
rate F'^'^'^ is defined as the inverse of the corresponding time scale. Assuming F'^"'^ ^ ljq, Ref. [20] derives a relation 
between decoherence and damping rates valid for arbitrary values of temperature: 



^dec 



icothf^ 
4 \2kBT 



( AZ 



V AZo 



(34) 



Formally, the factor coth (?ia;o/2fcBT') originates from the fluctuation-dissipation theorem, since decoherence may be 
described as a process of diffusion in phase space, which washes out the interference oscillations characterizing the 
coherence of the superposition state [22]. At high temperatures, this factor is approximately Ik^T j {J\io). Jointly with 
the high-temperature behavior of F, it results in dependence. This power law, also obtained in Ref. [23] in the 
case of a free particle, is visible in Fig. 4 (dashed line), where we plot F'^°'^/F^°^ = coth( ^l^) T/Fvac as a function of 
kBT/ (Tiujo). Both F and F'^'^'^ are already one order of magnitude larger than the vacuum values for ksT/ (hoja) ~ 0.4, 
showing that the thermal effect dominates even for moderate values of temperature. 

According to (34), decoherence is faster the larger is the separation AZ between the wave packets. For large 
separations, 'which way' information is more readily available since the two state components are better resolved, as 
compared to their widths in phase space. At zero temperature, zero-point fluctuations define the corresponding scale 
of length AZq. In the high-temperature limit, on the other hand, the additional factor 2kBT / (huio) in (34) yields 
pdcc^P ^ (AZ/Ar)^/2, the thermal de Broghe wavelength At = fi/ ^2mkBT replacing AZq as the reference of length 
scale. 

In conclusion, coherent superpositions of states very far apart in phase space are extremely unstable, since they 
decay into the corresponding statistical mixtures in a very short time. In the opposite extreme, coherent states are 
the most stable ones, in the sense that they generate the least entropy when interacting with the field [20]. They 
correspond to the 'pointer states', which play an important role in the classical limit of quantum mechanics [24]. Here 
we have discussed a particular physical mechanism of decoherence, which illustrates some of its general properties. 
The corresponding orders of magnitude are discussed in the end of the next section. 
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V. DISCUSSION AND FINAL REMARKS 



The temperature of the field defines the frequency scale ksT/h,. The motion is slow when the mechanical frequencies 
satisfy u) <C ksT/h. Hence the quasi static regime corresponds to the high-temperature limit. We summarize below 
the results obtained for this regime, now taking the time domain. 

For the scalar field, we have found 

The dissipative force results, in general, from input modes corresponding to regions i?2 and Rs in Fig. 1. R3 

corresponds to the input modes involved in the process of photon creation. At zero temperature, only R3 contributes, 
whereas in the quasi-static (high temperature) limit i?i and R2 provide the dominant contribution as shown in 
Appendix B. Thus, the dissipative component in (35) (first term in its r.-h.-s.) originates from scattering of thermal 
photons, rather than from creation of photons. It is of the form [25] —2Pt/c^ ^3(i)) where Pt, representing the total 
power incident on (both sides of) the mirror, is proportional to as predicted by the Stefan-Boltzmann law. This 
viscous force results from Doppler shifting the reflected thermal photons [9] . 

The dispersive component also appearing in (35) corresponds to a temperature-dependent, finite, and very small 
mass correction 

A™ = -fl(«2!A (36) 

27r n^c^ ^ ' 

For T = 300/r, we find Am/^ w — 1.5x 10~^^g/cm^. Am results from the contribution of input propagating modes that 

are scattered into evanescent waves (region i?4 in Fig. 1). It corrects the (experimentally known) zero-temperature 
mass of the mirror, which already contains tlie mass renormalization generated by vacuum fluctuations [14]- [16]. 
Even at zero temperature, the mass is modified when a second plate is present [26] [27] , or if its surface is corrugated 
[28]. In all three cases, the mass correction is a finite measurable quantity that depends on some control parameter: 
temperature, distance between the plates, or amplitude of corrugation. 

For the electromagnetic field, the viscous force is twice as large as in the scalar case, the two polarizations providing 
identical contributions: 

r^(t) = -^(^A5g(i). (37) 

Eq. (37) agrees with Ref. [20] , where the force is derived by considering the Doppler shift of the thermal photons. For 
T = 300K, the pressure is 4.1 x 10~^^N/m^ for a velocity of Im/s. 

For the dispersive force, on the other hand, the contributions from the two field polarizations cancel, as discussed 
in Sec. Ill, and the quasi-static force as given by (37) does not contain a term proportional to the second-order time 
derivative. Thus, in contrast to the scalar case, the mass correction vanishes in the electromagnetic case (it also 
vanishes for a scalar field in one dimension [9]). This difference between the two models is better understood by 
considering the thermodynamics of the field at temperature T in the presence of a static mirror. In Appendix C, we 
compute the Helmholtz free energy 

T =-kBT\a.Z, (38) 

where Z is the partition function, for both models. For the scalar field, the free energy in a volume V is = E^^^^T^^, 
where is the zero-point energy, and 

i,__^(fc^ C(3)(fcBr)^^ 

The first term in (39) is the free-space thermal free energy, featuring the characteristic T'* power law. For us, the 
interesting term is the second one. It is an effect of the Dirichlet boundary condition at the mirror's surface, which 
eliminates field modes propagating along directions parallel to the mirror. In fact, this term is not present in the 
electromagnetic case: 

yEM _ '^^ {ksTY 
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because of the contribution from TM-polarized grazing field modes (whereas the free-space term is twice the scalar 
result due to the equal contributions from the two orthogonal polarizations). 

From (39) and (40) we may also compute the entropy 5* = —dTj&T = —dJ^x/dT and the internal energy U = 

•D , rrD 
vac ' 1 



T + TS. For the scalar field, we find = + Z7^, with 



and MJt = Am c^. AUt represents the difference between the thermal energies in the presence of the static mirror 
and in free space. As discussed in Appendix D, where we present an alternative, 'local' derivation of U^, the energy 
density is decreased significantly (with respect to its free-space value) over a distance from the mirror of the order of 
hc/kBT. 

For the electromagnetic field, on the other hand, both !F^^ and are not modified by the mirror according to 
(40) (this may also be inferred from the result of Rcf. [29] for a two-plates configuration by taking the limit of large 
separation), and consistently the mass correction vanishes in this case. 

In the low temperature (or high mechanical frequency) limit, thermal photons provide a small correction of the 
vacuum dissipative force, of the order of {kgT/Tiuj)^ ^ 1, whereas the thermal dispersive correction is exponentially 
small, for both scalar and electromagnetic cases. We have also presented numerical results for arbitrary values of 
temperature, allowing us to discuss damping and decoherence of a mirror in a harmonic potential well in this general 
case. For both effects, thermal fiuctuations dominate over zero-point fluctuations for temperatures above 0.3?iLOo/kB, 
and the high temperature approximation provides accurate values for temperatures above hujo/ks- Although the effect 
of damping of energy is usually negligible, thermal radiation pressure efficiently destroys the coherence of a quantum 
superposition state. For T = 300K and A = Icm^, the decoherence time is l/r^^'^ = 1.3 x 10~^^s when the distance 
between the wave packet components is only AZ = Inm [30] . 
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APPENDIX A: DISPERSIVE THERMAL SUSCEPTIBILITY 

As discussed in Sec. II, the dispersive thermal susceptibility Rexr('-^) is given by the term proportional to n in 
Eq. (15), to be integrated over the region R4 in Fig. 1 (for the sake of clarity we assume w > in the derivation so 
as to be able to refer to Fig. 1): 

Rex?(w) = --^ y" J dwindk\\k\\^k'^^ - {uj + ui-.^Y ^ ui"^^ - /c^ n{-ujin) (Al) 
Eq. (Al) shows that RexT('^) < fo^' ^- We change to the variables 

^ _ n{k\\ - (Jin) _h{j^+^ 

ksT ' ksT ' ^ ' 



The Jacobian for this transformation yields 

|J|E 

Prom (Al) and (A3) we find 



d{u, v) 



d{k\\,Uin) 



ksT 



2 



6 2 — 1 



(A3) 



(A4) 



where O = Uu!/{kBT). We use (A4) to calculate Bjexri'^) numerically (see figure 2). Because of the exponential 
factor, the dominant contribution comes from values u ^ 1. In the high-temperature approximation, Q <C 1, we have 
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w < n <C u. Hence we may neglect and flv inside the root and v in the exponential in (A4), and replace the lower 
limit of integration over u by zero. We get 



Rex?M = -^(^) I dvV^T^) J^^ du-^^[l + 0{n')] (A5) 

The remaining integrals in (A5) give 27rC(3)f2^, yielding the result of (18). 

In the low-temperature limit, Q » 1, we may approximate the average photon number as follows: 

{exp [{u + v)/2] - 1}"^ « exp [-{u + v)/2] , 

because {u + v)/2> 12/2 in (A4). We change to new variables s and t with v/^ = and it/ri = 1 + and integrate 
over t by parts to find 

Rex?M«-^ (^M:y'j74g-n/2 dss^^l^e-""'"/^ p dt^ [tVl+^{l +f - s')]e-''''/'. (A6) 

The integral over t gives ^/2^^/^^{l — s^) by the method of steepest descent, the saddle point t = corresponding 
to the condition fc|| — cJin = lu. The resulting integral over s is also calculated by first integrating by parts, and then 
using the steepest descent method. The saddle point is at s = 0, corresponding to Win = — Thus, the dominant 
contribution in (Al) comes from the close neighborhood of = —Win = w/2, which corresponds to nearly grazing 
modes, and the final expression is given by (17). 

As discussed in Sec. Ill, the contribution of TE-polarized modes to the electromagnetic susceptibility turns out to 
be identical to The TM contribution is given by the r.-h.-s. of (23), its real T— dependent part yielding 

^ ' An^y n J Jo Jn - 1 ^uv{u - - v)' ^ ' 

It is clear from this equation that Re Xy^^^'' (w) is positive for any frequency w. Using the same methods employed in 
the scalar case, we derive from this equation the limiting cases corresponding to high and low temperatures. 
In the high- temperature limit, we derive from (A7) 

Re;,CTM) M /fc^ry i^-v)' r^^^n^0(^2)T (A8) 

^ ^ A7r^ \ n J Jo ^v{n - v) io et - 1 L ^ ^ ' 

The integrals above are readily calculated, resulting in (25). 

In the low temperature limit, we change to the variables s and t employed in the derivation of the scalar susceptibility, 
and derive from (A7) 

Rex<™'(.) » ^ {^)' 0^-" /; £ [(1 + *^)(»^ - 1/2) - m 

The integral over t gives y^27r/r2 (1 — s^)^/4 by the method of steepest descent (saddle point at t = 0). Then the 
integral over s is also computed with this method (saddle point at s = 0), yielding the result of (26). As in the scalar 
case, the dominant contribution comes from nearly grazing waves with frequencies close to —u/2. 



APPENDIX B: DISSIPATIVE THERMAL SUSCEPTIBILITY 

We calculate G^{t), defined in (19), from (15), taking the imaginary part of the term proportional to n. We change 
to the dimensionless variables il' = Tiu)i^/{kBT) and K = Tik\\/{kBT). For the contribution from region i?2 in Fig. 1, 
we also change f2' into —Q, — O', with Q. = Tiuj/{kBT) = 1/r. The joint contribution from J?i and i?2 is given by 

We split region into two sub-regions, corresponding to the intervals —uj < tOin < — w/2, and —lu/2 < Win < 0. We 
change SI' into Cl + Q' when integrating over the first sub-region, and into —fl' for the second sub-region. We find 
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Gf,){T) = 360r' dn'J^ dKK^in'^ - K^) [(O - - K^] (^^^73^ + _ J ■ (B2) 

In the high-temperature Hmit, ^ 1, we take e^^ — 1 w O and neglect il in the integrand in (Bl). We find 
<^(i+2)(t ^ 00) = 247rV. For the evaluation of Gf^y we take l/(e^^' - 1) 1/^]' [and likewise for l/{e^^-^' - 1)] in 
(B2) since SI' is bounded by 0/2. The final result is 0{t) and can be neglected. This completes the derivation of the 
high-temperatnrc limit of G^{t). 

In the low-temperature limit, O 1, we replace l/(e^+^ — 1) by e~^~^ , canceling the exponential factors in the 
numerator in (Bl). Actually, this is equivalent to neglect the contribution from R2. In fact, in this approximation 
only the close neighborhood of the origin in Fig. 1 contributes (as discussed in Appendix A, the dispersive component 
is exponentially small by a similar reason, since it results from region R4, which is far from the origin). We also 
take ^(O ± n')2 - K'^ w O in (Bl) and (B2), and in (B2) neglect l/(e"-"' - 1) and replace the upper bound of the 
integral over Cl' by infinity. In this approximation, the contributions from Ri and Rs are equal, and we find 

f°° Q'^ 

Gfi+2)(T-0) = Gf3)(T^0) = 12OTM dn'^^. (B3) 

Jo e — 1 

The remaining integral in (B3) gives 7r'*/15, completing the derivation of the low temperature limit of G^{t). 

In order to derive the dissipative susceptibility for the electromagnetic field, we need to consider the contribution 
of TM-polarized modes, which is given by the T— dependent imaginary part of the expression given by (23). The 
contributions from i?i and R2 give 

G ™Ut) = 360T5(e^^ - 1) / dSl' / dK , ^ ^ ^ ,^ ,a^a. (B4) 



whereas the contribution from R^ gives 

G^™\t) = 360t5 / dSl' / dK , " ^" ' " " " = ( — + 1 . (B5) 

In the high temperature limit, the contribution from Region R3 is 0(r) and hence negligible as in the scalar 
case, whereas the contribution from Ri and i?2 is derived from (B4) by neglecting O in the integrand and by taking 
e^^ - 1 « n. We find G^™\t ^ 00) = G(™)(^ ^ oo)[= G^{t ^ 00)]. 

For low temperatures, we proceed as in the discussion of the scalar case: the contributions from Ri and R3 are 
identical, and from R2 negligible. We find G'^™\t ^ 0) = 487r*T*. When added to G(™)(t ^ 0) = IStt^t^, we get 
the electromagnetic result as discussed in Sec. III. 



APPENDIX C: FIELD FREE ENERGY WITH A STATIC MIRROR 



We discuss both scalar (with Dirichlet boundary conditions) and electromagnetic field models. We first consider 
two parallel mirrors at rest separated by a distance L, and compute the free energy in the volume V = AL contained 
within the mirrors. Then, we take the limit X ^ 00 so as to identify the single-mirror effect. In this limit, we find a 
term proportional to V, corresponding to the free-space case, and, in the scalar case, a term proportional to A, which 
does not depend on L. The latter contains the independent single- mirror effects of the two mirrors, with only one side 
of each mirror taken into account. By symmetry, this is equal to the effect of one mirror when both sides arc taken 
into account. The free energy for two parallel plates has been calculated by several different methods [31]. Here we 
follow the approach of Ref. [32] . 

The partition function is defined as 

Z = tvexp[-Ho/{kBT)], (CI) 

with Ho = ?kv\{a\ax + 1/2). A = {n, ky} represents a set of parameters defining a given field mode (which also 
contains an index for polarization in the electromagnetic case), with ux = ^/{mr/L)'^ + k?,. Eq. (CI) yields 



1 _ e-fta,x/(fcBT) ' (<^2) 
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which combined with (38) leads to 



^ = ^ ^ + fcsT^ In (l - e-'^--/('=-^)) . (C3) 



The first term in (C3) represents the zero point energy fyac- We want to compute the second term, which contains 
the temperature dependence: 

oo 

TT = Y,g{n)a{n), (C4) 

n=0 

where 



ayn) = — / a«;||A;||ln 



1 — exp 



(C5) 



with g{n) = 1 — (5no and g{n) = 2 — (5„o in the scalar and electromagnetic models, respectively. Hence in the former case 
grazing modes, corresponding to n = 0, are ruled out, whereas in the electromagnetic model TM polarized grazing 
modes provide a non vanishing contribution to the free energy. 
In order to compute Tt, we employ Poisson sum formula 

^ oc oo 

-a(0) + ^o(n) =Ao + 2^ReAj, (C6) 

71=1 j=l 

where 

/■oo 

Aj = / rfna(n)e*2-'''^". (C7) 
Jo 

We change the variable of integration from k^^ to ^ = Ti^ [n-K / LY + kyiksT), and integrate over n to find 

^°-"90^^' ^^^^ 



and for non vanishing values of j, 



27r2jTi2 Jo 2i 



^i = ^T^ d^Cln(l-e-€)^— ^, (C9) 



where = 2jLkBT/U. By integrating the r.-h.-s. of (C9) by parts, we find 

AiksTf j 



TipA ■ — 

' ~ Sir^phL 1 2A 



^ [ttAj coth(7rAj) - 1] + Re ^ '^T^'"'^'^ } ' 



We calculate the r.-h.-s. of (CIO) in the limit A^ 3> 1. By successive integrations by parts, one may show that the 
last term within brackets in (CIO) yields — 1/(2A|) + 0{l/A^) and may be neglected in (CIO): 

ReA, = -^^^[1 + 0{1/Aj)]. (Cll) 

Hence in the limit L — > oo only contributes in the Poisson sum formula Eq. (C6). 

Finally, when computing the free energy for the scalar Dirichlct field, we need to subtract a(0)/2 from the r.-h.-s. 
of (C6), since (?(0) = in (C4) in this case. We calculate a(0) from (C5): 

a(0) = -^^li^. (C12) 

Combining (C4), (C6), (C8) and (C12), we derive the free energy for the Dirichlet field and a single static plate as 

given by (39). For the electromagnetic field, on the other hand, wc may employ Poisson formula directly, since in this 
case the sum over n in (C4) has already the form of the l.-h.-s. of (C6). Then we find !F^^'^ = 2Ao [with Aq given by 
(C8)], in agreement with (40). 
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APPENDIX D: THERMAL ENERGY DENSITY FOR THE DIRICHLET SCALAR FIELD WITH A 

STATIC MIRROR 



In this appendix, we compute the energy density 



(Dl) 



for a scalar field at temperature T, and with Dirichlet boundary conditions at z = 0. Since we consider the static case, 
we replace </> by (/)°, and from its normal mode expansion as given by Eq. (7) derive (with k = y'fc| + fcf ) 



w{z) = 



dk. 



dkn 



{kP' + fcjj^ siv?{kzz) + kl cos^{kzz) 



(D2) 



As in the derivation of the susceptibility xi^)} the energy density is naturally split into two contributions, one from 
vacuum fluctuations [the '1/2' inside brackets in (D2)], the other from thermal fluctuations, which corresponds to the 
factor n(fc) in (D2). Accordingly, wc write w = wt + Wvac- Here wc arc only interested in the thermal contribution 
wt, which is obtained from (D2) after changing the variable of integration from fcii to k : 



wt{z) = — ^3 + AwT. 



(D3) 



The first term in (D3) represents the free-space energy density for a scalar field; the effect of the boundary at z = 
is contained in Awt ■ 



Awriz) = ^ 



J__7rfcBr f 2TrkBTz \ ^ ( kpTV 1 

2^4 2 nz^ V ^ ) ) sinh^ {^^''f") 



(D4) 



Awt{z) is a negative-defined, increasing function of z that goes to zero as —kTj (167r2;^) for z ^ h/{kBT). Hence the 
mirror reduces the thermal energy density, an effect stronger near the mirror. Awt is finite at ^; = 0, and vanishes at 
T = as expected from its definition. 

The total modification of the internal energy is given by the volume integral of Awt{z), taking into account both 
sides of the mirror: 



AUt = 2A I dz Awt{z) = 



f 

O — < 



dZ 



2 coth.E 
'2A 



Z'^ sinh^ Z 



(D5) 



The integral above may be computed with the method of residues, by taking a semi-circular path (radius oo) in 
the complex plane of Z. The poles of the integrand lie along the imaginary axis, at the positions Z = inir, with n 
integer, n^O. We find 



ATI - 1 



Z-/ r) 3 



(D6) 



The series appearing in (D6) is equal to C(3)- Then, comparing (36) with (D6), we find AUt = Amc? as discussed in 
Sec. V. 
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